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ABSTRACT 

A few particles of presolar AI 2 O 3 grains with sizes above 0.5 are believed to have been produced 
in the ejecta of core-collapse supernovae (SNe). In order to clarify the formation condition of such 
large AI 2 O 3 grains, we investigate the condensation of AI 2 O 3 grains for wide ranges of the gas density 
and cooling rate. We first show that the average radius and condensation efficiency of newly formed 
AI 2 O 3 grains are successfully described by a non-dimensional quantity Aon defined as the ratio of 
the timescale with which the supersaturation ratio increases to the collision timescale of reactant gas 
species at dust formation. Then, we find that the formation of submicron-sized AI 2 O 3 grains requires 
at least ten times higher gas densities than those presented by one-dimensional SN models. This 
indicates that presolar AI 2 O 3 grains identified as a SN origin might be formed in dense gas clumps, 
allowing us to propose that the measured sizes of presolar grains can be a powerful tool to constrain 
the physical conditions in which they formed. We also briefly discuss the survival of newly formed 
AI 2 O 3 grains against the destruction in the shocked gas within the SN remnants. 

Subject headings: dust, extinction - Galaxy: evolution - ISM: supernova remnants - meteorites, 
meteors, meteoroids - stars: massive - supernovae: general 


1. INTRODUCTION 

The role of core-collapse supernovae (SNe) as dust pro¬ 
ducers is a fundamental issue for understanding the evo¬ 
lution history of dust in the universe. Recent far-infrared 
observations of young supernova remnants (SNRs) re¬ 
vealed that 0 . 1-1 Mq of dust had formed in the metal- 
rich inner ejecta (Sibthorpe et al. 2010; Barlow et al. 
2010; Gomez et al. 2012; Matsuura et al. 2011, 2015; 
Indebetouw et al. 2014). However, it remains to be 
clarified what fraction of these newly formed grains can 
survive the destruction in the shock-heated gas within 
the SNRs and be injected into the interstellar medium 
(ISM). The destruction efficiency of dust grains depends 
on their chemical compositions and size distributions 
(e.g., Nozawa et al. 2007), which are determined by the 
density and temperature evolution of the gas out of which 
they form (Nozawa & Kozasa 2013) as well as the degree 
of mixing and dumpiness of the ejecta (Nozawa et al. 
2003; Sarangi & Cherchneff 2015). Therefore, the phys¬ 
ical condition and structure of the SN ejecta must be 
fully appreciated through various approaches to unravel 
the properties of dust that is formed in the ejecta and is 
finally ejected to the ISM. 

Presolar grains, which are identified in meteorites due 
to their highly anomalous isotopic compositions, are in¬ 
valuable fossils that enable us to directly observe the 
detailed chemical compositions and sizes of stellar dust 
(Clayton & Nittler 2004 and references therein). Their 
isotopic signatures give clues about nucleosynthesis in 
stars and mixing of elements in the SN ejecta. In addi- 
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tion, the measured sizes of presolar grains could offer key 
information on the physical conditions at their formation 
sites. 

Among the presolar grains which are considered to 
have originated in SNe, AI 2 O 3 grains are of great impor¬ 
tance because AI 2 O 3 is believed to be one of the major 
dust components in Cassiopeia A (Gas A) SNR (Douvion 
et al. 2001; Rho et al. 2008). Furthermore, most of the 
dust formation calculations have predicted the formation 
of AI 2 O 3 grains in the ejecta of SNe as the first conden¬ 
sate among oxide grains (Kozasa et al. 1989, 1991; Todini 
& Ferrara 2001; Nozawa et al. 2003, 2008, 2010; Bianchi 
& Schneider 2007; Kozasa et al. 2009; Sarangi & Cherch¬ 
neff 2015). However, the calculated sizes of AI 2 O 3 grains 
are below ^0.03 ^m, which is much smaller than the 
measured sizes of presolar AI 2 O 3 grains with 0.5 ^m to 
1.5 p,m in diameter (Nittler et al. 1998; Choi et al. 1998), 
and such small AI 2 O 3 grains are found to be almost com¬ 
pletely destroyed in the hot gas within the SNRs before 
being ejected to the ISM (Nozawa et al. 2007; Silvia et 
al. 2010, 2012). This seems to contradict the fact that 
we are observing the AI 2 O 3 grains of a SN origin in hand 
on the Earth. 

One of the main reasons that only small AI 2 O 3 grains 
are produced in the simulations is the low number density 
of Al atoms in the ejecta, led by the relatively homoge¬ 
neous ejecta of spherically symmetric SN models (here¬ 
after referred to as 1-D SN models). In reality, the SN 
ejecta should be much more inhomogeneous and compli¬ 
cated, as is suggested from multi-dimensional hydrody¬ 
namic simulations (e.g., Kifonidis et al. 2003; Hammer 
et al. 2010; Joggerst et al. 2010) and various observa¬ 
tions of SN 1987A (e.g., Kj*r et al. 2010; Larsson et al. 
2013). This implies that the formation of dust grains 
would proceed in the gases with a variety of densities, 
and that large grains could be formed in high-density 
clumps. 
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In this Letter, we investigate the formation of AI 2 O 3 
grains for wide ranges of density and cooling rate of 
gas, to explore the formation condition of presolar AI 2 O 3 
grains as large as measured in meteorites. We show that 
submicron-sized AI 2 O 3 grains can be produced only in 
the gas with more than ten times higher densities than 
those predicted by 1-D SN models, suggestive of the pres¬ 
ence of dense gas clumps in the ejecta. We also examine 
the survival of AI 2 O 3 grains against the destruction pro¬ 
cess in SNRs. We propose that the comparison between 
the calculated sizes and the measured sizes of presolar 
grains can be a novel and valuable approach that gives 
insight into the physical conditions and structure of the 
SN ejecta. 

2. MODEL OF DUST FORMATION 

The formation of AI 2 O 3 grains is calculated by apply¬ 
ing the formula of non-steady-state dust formation in 
Nozawa & Kozasa (2013). In this formula, the forma¬ 
tion of small clusters and the growth of grains are self- 
consistently followed under the consideration that the ki¬ 
netics of dust formation process is controlled by collisions 
of key species, defined as the gas species that have the 
lowest collisional frequency among the reactants. The 
formula leads us to derive the size distribution and con¬ 
densation efficiency of newly formed grains, given chem¬ 
ical reactions for the formation of clusters, abundances 
of the relevant gas species, and time evolutions of gas 
density and temperature. The detailed prescription of 
the calculations of non-steady dust formation is given in 
Nozawa & Kozasa (2013). 

In the ejecta of SNe, the most likely formation site of 
AI 2 O 3 grains is the 0-rich layer, where Al atoms, as well 
as O atoms, abundantly exist. We consider as a chemical 
reaction at cluster formation 2A1 -|- 30 ^ AI 2 O 3 (Kozasa 
et al. 1989). The initial number ratio of Al to O atoms 
is taken as cai.o/co,o = 1 / 200 , where cai,o and co,o are, 
respectively, the number densities of Al and O atoms at 
a given initial time t = to, so that the key species are Al 
atoms. This abundance ratio approximately corresponds 
to that in the Al-rich region of solar-metallicity SNe (see, 
e.g., Kozasa et al. 2009; Nozawa et al. 2010). Note that 
the results of calculations are little affected by the Al/0 
ratio as long as cai,o/co,o 1 - 

The ejecta of SNe freely expands since ^1 day post¬ 
explosion, and the gas density is inversely proportional 
to the cube of time t. Thus, the number density of a gas 
species Ci{t) (where i is Al or O), without the depletion 
of the gas-phase atoms due to the formation of clusters 
and grains, is given by 

Ci(i) = Ci,o ■ (1) 

As in Nozawa & Kozasa (2013), the gas temperature T{t) 
is assumed to decrease as 

-3(7-1) 

T«) = r.(-) . (2) 

where Tq is the gas temperature at tp, and 7 is a free 
parameter that prescribes the cooling rate. 

As the gas cools down, the formation of dust from gas 
can be realized in a supersaturated state (S' > 1); S' is 
the supersaturation ratio defined as In S = —Ag/kT with 


TABLE 1 

Numerical Constants Used for Dust Formation 
Calculations 


A/IO"* K 

B 

ao^ 




(A) 

(erg cm“^) 

18.4788 

45.3542 

1.718 

690 


^hypothetical radius of the condensate per key molecule. 
I’surface tension of bulk grains (Overbury et al. 1975). 


k being the Boltzmann constant and Ag the change of 
the chemical potential per key species for the formation 
of bulk condensate from the reactants, as is fomulated in 
Equation (30) in Nozawa & Kozasa (2013). Here we take 
to as a time at which S = 1 and determine Tq for given 
CAi.o and co,o from the equation 


lnS = - 


kTn 


In 


CAl.ofc'To 


Ps 


■In 


cofikTo 

Ps 


= 0,(3) 


where Ag is the change of the chemical potential at the 
standrad pressure Ps and is approximated as Ag/kT = 
—A/T+B with the numerical values A and B taken from 
Nozawa et al. (2003) (see Table 1). 

In this study, we suppose the clusters containing more 
than 100 Al atoms as bulk grains. This corresponds to 
the minimum grain radius of a* = 8.0 A (AI 2 O 3 grains 
are assumed to be spherical). The sticking probability of 
gas species is assumed to be unity for any sizes of clus¬ 
ters and grains. In the calculations, we take to = 300 
days, so the free parameters are cai,o and 7 , for which 
we consider the ranges of cai,o = 10 ^- 10 ^^ cm“^ and 
7 = 1.1-1.7. The calculations are performed until the 
gas density becomes so low that grain growth is negligi¬ 
ble. In what follows, we mainly examine the resultant 
behavior of the final average radius Cave.oo and conden¬ 
sation efficiency /con ,00 that are obtained at the end of 
the calculations. The condensation efficiency is defined 
as the fraction of Al atoms locked up in AI 2 O 3 grains. 


3. RESULTS OF DUST FORMATION CALCULATIONS 

Figure 1 shows the size distributions of newly formed 
AI 2 O 3 grains calculated for 7 = 1.25, adopting cai,o = 
10^, 10®, and 10® cm“®. For these cases, all Al atoms are 
finally locked up in AI 2 O 3 grains (that is, /con,ex, = !)■ As 
seen from the figure, the size distribution is lognormal¬ 
like for any of cai,o considered here, with a narrower 
width for a higher cai,o- More importantly, the average 
radius increases with increasing cai,o: aave,oo = 0.0067, 
0.047, and 0.36 fim for cai,o = 10^, 10®, and 10® cm“®, 
respectively. This is because a higher gas density leads 
to more efficient growth of grains. The results in Figure 
1 point out that, for 7 = 1.25, the number density of Al 
atoms at dust formation must be higher than ~5 x 10 ® 
cm“® in order that AI 2 O 3 grains with the radii larger 
than 0.25 pm (^0.5 pm in diameter) can be formed. 

Nozawa & Kozasa (2013) demonstrated that the for¬ 
mation process of dust grains is described in terms of the 
timescales of two physical quantities: the timescale with 
which the supersaturation ratio S increases Tsat and the 
collision timescale of key species Tcqu. They found that, 
for C and MgSiOR grains, the average radius and con¬ 
densation efficiency are universally scaled by one non- 
dimensional quantity Aon = Tsat(ton)/Tcon(ton), where 
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Fig. 1. — Final size distribution spectrum of newly formed AI 2 O 3 
grains calculated for 7 = 1.25 and to = 300 days. The size distribu¬ 
tions are plotted as a^n{a) so as to represent the mass distribution 
per logarithmic grain radius, where n(a) is the arbitrarily normal¬ 
ized number density of grains with radii between a and a + da. 
Three cases are considered for the number density of A1 atoms at 
t = to; cai,o = (dotted), 10 ® (solid), and 10 ® cm“® (dashed). 

ton is the onset time of dust formation (ton ^ to) and 
is taken as a time at which the condensation efficiency 
reaches 10“^°. According to their study, it would be in¬ 
teresting to see if such a scaling relation holds for AI 2 O 3 
grains. 

Figure 2 depicts the average grain radii Oave.oo and con¬ 
densation efficiencies /con ,00 obtained from the calcula¬ 
tions with a variety of cai.o for each of 7 = 1.1, 1.3, 1.5, 
and 1.7, as a function of Aon- In the present study, Aon 
is approximately written as (Nozawa & Kozasa 2013) 


/ ton \ 

CAl(^on) 

’ T{ton) ' 

\300 days/ 

10 ® cm“® 

2000 K 


Figure 2 clearly shows that the results for different 7 are 
well overplotted, indicating that both Cave.oo and /con.00 
are uniquely determined by Aon even for AI 2 O 3 grains. 
As is the cases of C and MgSi 03 grains, the formation of 
AI 2 O 3 grains can be realized at Aon ^ 1 “ 2 , and /con .00 = 
1 at Aon > 20 . In addition, aave.oo becomes large as Aon 
increases, which means that the final average radius is 
larger for a higher gas density and/or a slower gas cooling 
because Aon is roughly proportional to the product of gas 
density and cooling timescale that is reflected by ton- 

In Figure 2, we also plot the ranges of Aon expected 
in the ejecta of Type ITP and 11b SNe, referring to 1-D 
SN models used in Kozasa et al. (2009) and Nozawa et 
al. (2010), respectively. For a Type ITP SN with the 
massive hydrogen envelope, the number density of Al 
atoms in the Al-rich region is estimated to be CAi(ton) — 
(0.2-8) xlO”^ cm“^ at ton = 300 days. For 7 ~ 1.25 and 
T{ton) — 2000 K, this corresponds to Aon — 70-3000, for 
which Oave.oo — 0.002-0.03 pm. On the other hand, for 
a Type 11b SN with the small-mass envelope, CAi(ton) — 
(0.3-5) xlO^ cm“^ at ton = 300 days, resulting in Aon — 
0.1-2. Hence, AI 2 O 3 grains are not expected to form 
in the Type lib SN. Nozawa et al. (2010) reported the 
formation of AI 2 O 3 grains in the Type Ilb SN but their 
average radii are less than A, which is regarded as 
small clusters in this study. 

Our calculations show that, in order to produce AI 2 O 3 


^on ”^sat(^on^/^coll(^on) 

Fig. 2.— Dependence of the final average radii aave.cxj and con¬ 
densation efiiciencies /con,00 of newly formed AI 2 O 3 grains on Aon- 
The results for four different 7 (7 = 1.1, 1.3, 1.5, and 1.7) are shown 
in each color but they are plotted as almost the same curve. The 
hatched regions depict the expected ranges of Aon for the forma¬ 
tion of AI 2 O 3 grains in the Al-rich region, referring to 1-D models 
of a Type II-P SN (cyan; Kozasa et al. 2009) and a Type lib SN 
(Nozawa et al. 2010). The solid vertical line denotes the minimum 
value of Aon (= 3 X 10^) necessary for explaining the measured 
sizes (radius of >0.25 fim) of presolar AI 2 O 3 grains. 

grains larger than 0.25 pm (diameter of >0.5 pm) as 
measured for presolar grains, Aon should be higher than 
3 X 10"^. Given that 7 and T(ton) do not change largely, 
such a high Aon could be achieved by considering the 
gas densities that are more than one order of magnitude 
higher than those presented by 1-D SN models (namely, 
CAi(ton) ^ 7 X 10® cm“® at ton = 300 days, see Equation 
(4)). This strongly suggests that the discovered presolar 
AI 2 O 3 grains were formed in dense clumps within the 
ejecta. Equation (4) tells us that the formation of large 
AI 2 O 3 grains may also be possible if the formation time 
of dust is later than ton = 3000 days. However, it is 
too hard to keep the gas density as high as cai — 10 ® 
cm“® in such late epochs because the gas density rapidly 
decreases with time (see Equation (2)). 

4. CONCLUSIONS AND DISCUSSION 

We have systematically investigated the formation of 
AI 2 O 3 grains, adopting a wide variety of gas densities 
and cooling rates. We show that, as in the cases of 
G and MgSiOs grains, the average radius and conden¬ 
sation efficiency of AI 2 O 3 grains are nicely scaled by a 
non-dimensional quantity Aon defined as the ratio be¬ 
tween the timescale of the supersaturation ratio and the 
collision timescale of key species at dust formation. We 
also find that large AI 2 O 3 grains with radii of >0.25 pm 
can be formed only in dense gas regions which have more 
than ten times higher densities than those estimated from 
1-D SN models. This points out the presence of dense 
clumps in the ejecta of core-collapse SNe. 

The formation of dust in dense clumps was deduced 
from various early-phase (<1000 days after explosion) 
observations of SN 1987A (e.g., Lucy et al. 1991; Meikle 
et al. 1993; Golgan et al. 1994), which can be classified as 
a Type ITP SN. The recent radiative transfer models of 
dust emission and absoption also suggested the necessity 
of optically thick clumps to account for the evolution of 
optical to infrared emission from SN 1987A over ~600- 
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9000 days (Ercolano et al. 2007; Wesson et al. 2015; Dwek 
& Arendt 2015; Bevan & Barlow 2015). The density 
contrast between the clumps and interclumps considered 
by these works is ~ 10 - 100 , which is in good agreement 
with the density enhancement needed for the formation 
of large AI 2 O 3 grains (i.e., more than ten time that of 
1-D SN models). This allows us to conclude that the 
measured size of presolar grains is an independent and 
useful probe to constrain the clump density in the ejecta. 

Our calculations show that AI 2 O 3 grains cannot form 
in Type Ilb SNe as a result of the too low ejecta density 
in the 1-D SN model. However, the analyses of infrared 
emission spectra of Cas A, which was identified as Type 
lib through the light echo (Krause et al. 2008), have sug¬ 
gested the presense of newly formed AI 2 O 3 grains (Dou- 
vion et al. 2001; Rho et al. 2008). This contradiction can 
also be resolved by considering the dense gas clumps in 
the ejecta; from Figure 2, we can estimate that the den¬ 
sity of gas clumps required for the formation of AI 2 O 3 
grains (that is, Aon ^ 2 ) is more than ten times that of 
the 1-D Type lib SN model. We also note that the re¬ 
cent observation of near-infrared extinction implies the 
existence of large (>0.I ^m) Si grains in Cas A, provid¬ 
ing another indication of dust formation in dense clumps 
(Lee et al. 2015). 

Next we discuss the survival of AI 2 O 3 grains formed 
in the ejecta against the destruction by SN shocks on 
the basis of a dust evolution model in SNRs (Nozawa 
et al. 2007). This model assumes a spherically symme¬ 
try shock but it has been shown that the destruction 
efficiency of dust is not greatly different from those by 
multi-dimensional hydrodynamic simulations (Silvia et 
al. 2010, 2012). As the initial condition of the ejecta, we 
here consider the core-envelope structure of 3 -I- 12 Mq 
(implicitly assuming Type ITP SNe). The other input 
parameters are the explosion energy and the gas den¬ 
sity in the ISM, for which we take reasonable values of 
1.8 X 10^^ erg and I.O cm“^, respectively. 

The result of the calculation is given in Figure 3, which 
shows that small grains with the initial radii of Oini ^ 
0.01 ^m are completely destroyed in the shocked gas. 
Grains with 0.02 fim < Oini < 0.1 fim are eroded in 
the hot gas and are finally destroyed as soon as they 
encounter the cool dense shell that is formed behind the 
forward shock after 5 x 10"^ years. On the other hand, 
AI 2 O 3 grains with Oini ^ 0.2 /rm can be ejected from 
SNe without reducing their size significantly. This simple 
calculation illustrates that, once large AI 2 O 3 grains are 
produced in dense clumps, they are likely to survive the 
destruction in SNRs and to be easily transported to the 
ISM0Thus, we suggest that the formation of submicron¬ 
sized AI 2 O 3 grains in dense clumps is also indispensable 
in order that newly formed AI 2 O 3 grains can endure the 
destruction by shocks and be incorporated into nearby 
molecular clouds and protoplanetary disks. 

The identification of AI 2 O 3 grains as a SN origin comes 
from their relatively high ratio (so-called Group 

4 grains, Nittler et al. 1997; Choi et al. 1998) or large en¬ 
richment of (represented by Grain T84, Nittler et al. 

® In this calculation, we do not include the destruction of dust 
by non-thermal sputtering resulting from the high-velocity motion 
of dust after being injected into the ISM, in order to focus on the 
survivability of dust in the SNR. 



time (yr) 

Fig. 3. — Time evolutions of the positions (upper) and the ratios 
of radii to the initial ones (lower) of AI 2 O 3 grains in a SNR. The 
results are plotted for the initial grain radii of 0.01, 0.02, 0.05, 0.1, 
and 0.25 fim. In the upper panel, the thick solid curves depict the 
trajectories of the forward shock (FS) and reverse shock (RS). 

1998). Since we consider the formation of AI 2 O 3 grains 
in the 0 -rich layer where is rich, they may be cat¬ 
egorized as T84-like grains (or Group 3 grains with the 
moderate enhancement of ^®0). On the other hand, to 
explain the ^®0-enriched composition of Group 4 grains, 
the extensive mixing between different layers in the ejecta 
is invoked with the relatively 0 -rich hydrogen enve¬ 
lope being the dominant component (Nittler et al. 2008). 
However, even if the large-scale mixing takes place, the 
molecular diffusion lengths are much smaller than the 
typical size of gas clumps so that the microscopic mixing 
of elements may be very ineffective (Deneault et al. 2003). 
Therefore, the origin of oxygen isotopic composition of 
the Group 4 AI 2 O 3 grains, as well as their formation pro¬ 
cess, is still a challenging problem. 

We conclude that dense gas clumps are necessary for 
the formation of submicron-sized AI 2 O 3 grains as discov¬ 
ered in meteorites. The presence of such dense clumps 
in the 0 -rich layer may also cause silicate grains to be 
formed with very large radii, compared to those based on 
I-D SN models. Given that the radii of silicate grains are 
generally larger than those of AI 2 O 3 grains by a factor 
of about ten (e.g., Nozawa et al. 2003), the discovery of 
^®0-enriched micron-sized (~I-I0 /iin) silicate grains in 
meteorites will serve as further evidence for dense clumps 
in the SN ejecta. 
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